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Abstract 

The flavor dependence of jet quenching is a powerful handle to discriminate between models of 
parton energy loss in heavy-ion collisions. We demonstrate the capacity of CMS to identify jets 
initiated by bottom quarks using displaced vertices reconstructed in the silicon tracking system. 
The b-jet to inclusive jet ratio is measured in PbPb collisions and compared to pp collisions and 
simulations at the same center-of-mass energy. 



Identification of reconstructed jets from fragmentation of b quarks has so far not been per- 
formed in heavy-ion collisions. However, recent CMS data demonstrate that non-prompt J/ift 
mesons, i.e., those coming from decay of B hadrons, are indeed suppressed in PbPb collisions 
with respect to the pp expectation JT]. This indirect measurement, which samples typical ]/if/ pj 
of values on the order of 10 GeV/c, provides a strong motivation to perform a more direct mea- 
surement using fully reconstructed jets. Such a measurement would enable a direct comparison 
of b quark energy loss to that of inclusive jets at much larger values of jet pj, where the flavor 
dependence of parton energy loss can be probed in detail. 

CMS is described in detail in |2|. Jets formed from heavy flavor quark fragmentation can be 
tagged by the presence of displaced vertices, either by direct reconstruction of these vertices or 
by the impact parameter (i.e., the distance of closest approach to the primary vertex) of tracks 
originating from these vertices Q. Information from these tracks and vertices are typically 
combined into a quantity which optimizes their discrimination between heavy and light flavor 
jets. In this analysis, we use a discriminator to tag b-jets which is based on the flight distance 
of the reconstructed secondary vertex (SV) with respect to the primary vertex of the interaction. 
The efficiency of this SV tagging is then evaluated both directly from simulation, and with a 
data-driven method using a weakly correlated tagger based on the impact parameter of tracks 
associated to the jets. 

The performance of lifetime-based tagging relies on the high efficiency and low fake rate of 
reconstruction of charge particle tracks from displaced vertices. The standard heavy-ion tracking 
algorithm ||4) in CMS is largely restricted to the reconstruction of charged particles from the pri- 
mary vertex. To enhance the efficiency of tracks from secondary vertices, additional track recon- 
struction is performed. Secondary track reconstruction in central heavy-ion events is extremely 
resource intensive, due to the large number of possible hit combinations when incorporating hits 
over a sufficiently large region to find such displaced tracks. To mitigate the number of possible 
track candidates, the additional track reconstruction uses reconstructed jets as seeds and limits 
the search window for tracker hits to a region defined around the jet axis. 

Identification of b-jets is achieved using a discriminating variable, denoted Simple Secondary 
Vertex High Efficiency (SSVHE) Q, which is based on the flight distance of reconstructed S Vs. 
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Detailed studies of quantities related to SV reconstruction show that it is well modeled by PbPb 
simulations, lending confidence to the use of this discriminator to identify b-jets. Another dis- 
criminating variable, called Jet Probability (JP), is also used in this analysis. The JP algorithm 
orders jet-associated tracks based on their impact parameter significance (impact parameter di- 
vided by its uncertainty) and calculates the likelihood that they come from the same primary 
vertex J3J. The JP algorithm provides a measure of discrimination for nearly all b-jets, even in 
the case when no SV is reconstructed. This property of the JP tagger is exploited to obtain a 
data-driven estimate of the SSVHE tagging efficiency 0. This is done by evaluating the fraction 
of jets from bottom fragmentation with and without the SSVHE tagging requirement, fitting the 
JP discriminator distribution with the same template fitting method described below. 

The performance of the discriminators is typically benchmarked by plotting the b-jet tagging 
efficiency vs. the light (udsg) jet mis-tag efficiency, such that the resulting curves are indepen- 
dent of the underlying b-jet fraction. Fig. [T] shows the b-jet tagging efficiency vs. the light jet 
mis-tag efficiency from Pythia and Pythia+Hydjet. No centrality selection is applied to the 
Pythia+Hydjet sample, but its centrality distribution is weighted to match the distribution in 
data after the jet selection. The two taggers used in this analysis, SSVHE and JP, are shown. 
The performance of the b-tagging degrades somewhat with the increased multiplicity of PbPb 
collisions, giving roughly a factor of 3 poorer rejection of light jets for a b-jet efficiency of 50%, 
relative to Pythia alone. Despite the reduced performance, one is still able to achieve roughly 
a factor of 100 rejection of light jets for a b-jet efficiency around 50%. The charm rejection for 
this b-jet efficiency is about a factor of 10 (not shown). 

The ratio of b-jet to inclusive jets (the "b-jet fraction") is extracted from the data using the 
SV mass. The SV mass, which is simply the invariant mass of the charged tracks used in the SV 
reconstruction, provides good separation between the charm and bottom contributions. Unbinned 
maximum likelihood fits are then performed where the b-jet and non-b-jet shapes are fixed by the 
simulation, but their normalizations are allowed to float. Figure[2]shows S V mass distributions in 
data and corresponding simulation templates for two of the jet pj bins used in this analysis. The 
ratio of the charm to light normalizations are fixed in the template fits according to their relative 




Figure 1: The b-jet tagging efficiency vs. the light jet mis-tag efficiency for simulated pp events from Pythia (left) and 
simulated PbPb events from Pythia embedded in Hydjet (right) for the SSVHE and JP discriminators. The red cross 
marks the working point of the SSVHE discriminator used in this analysis. 
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Figure 2: Template fits to the SV mass distributions in PbPb collisions, after tagging with the SSVHE discriminator for 
jets of 80 < p T < GeV/c (left) and 120 < p T < 150 GeV/c (right). 



contributions in the input simulation. 

The SSVHE tagging efficiency is shown in the left panel of Fig. [3] The result obtained using 
JP as a reference tagger is compared to the result obtained directly from simulation. The b-jet 
purity is the fraction of jets from bottom quarks in the SS VHE-tagged sample. The right panel 
of Fig. [3] shows the b-jet purity of the tagged sample as a function of the jet pj extracted from 
the template fits in Fig. [2] These purity values are close to those from the input simulation. 

The b-jet to inclusive jet ratio is calculated from (A^ t ^ sed /A^ e ts)(^ > /e), where N^ ed and iVj ets 
are the number of tagged jets and the total number of jets counted from the data, respectively, 
and P and e correspond to the purity and efficiency plotted in the left and right panels of Fig. [3] 
respectively. The left panel of Fig. |4] shows this ratio as function of jet pj in 0-100% PbPb cen- 
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Figure 3 : Left: The tagging efficiency of the SSVHE discriminator from simulation and from the reference tagger method. 
Right: The b-jet purity extracted from template fits to the SV mass distributions, compared to the input simulation. 
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Figure 4: The b-jet to inclusive jet ratio in 0-100% PbPb collisions (left) and pp collisions (right) as a function of jet pi 
compared to Pythia embedded in Hydjet (PbPb) and Pythia (pp). 



trality collisions. The b-jet fraction is around 2.9-3.5% with no significant pi dependence, with 
an absolute uncertainty in the range of 0.6-1.1%, increasing with jet p T . The values predicted by 
Pythia+Hydjet, which decrease from about 3% to 2.5% as a function of pj, are consistent with 
the data. The b-jet fraction in pp collisions as a function of jet pj is shown in the right panel of 
Fig. [4] The values predicted by Pythia are consistent with the data. 

There are several sources of systematic uncertainty which affect the b-jet fraction. The un- 
certainty on the tagging efficiency is evaluated from the difference between the simulation and 
data-driven efficiencies. A second category of uncertainty is composed of effects which influ- 
ence the b-tagging purity. One such uncertainty is the relative contribution from charm and light 
jets, which changes the shape of the non-b-jet template. The relative normalization is fixed in 
our template fits, and then alternately allowed to float to obtain a systematic uncertainty. To de- 
termine the uncertainty on the modeling of the non-b-jet template shapes by the simulation, we 
compared to results obtained with a data-driven template formed by requiring an anti-b-tag on the 
JP discriminator. An uncertainty on the stability of the template shapes was estimated by varying 
the working point of the SSVHE tagger. Other sources of systematic uncertainty considered are 
the b-jet and inclusive jet energy scales and the modeling of b-jet production in Pythia. 

To within the fairly sizable uncertainties, the data disfavor an extreme scenario in which b-jets 
suffer no energy loss in PbPb collisions. A more detailed understanding of the jet p T dependence 
of the b-jet to inclusive jet ratio and, more generally, the parton mass and flavor dependence of 
energy loss, would require a reduction of the statistical and systematic uncertainties. 
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